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Abstract. We present simulations of the interferometric visibilities of Herbig Ae star disks. We investigate whether 
interferometric measurements in the 10 /im atmospheric window are sensitive to the presence of an increased scale 
height at the inner disk edge, predicted by recent models. Furthermore, we investigate whether such measurements 
can discriminate between disks with a "flaring" geometry and disks with a "flat" geometry. We show that both 
these questions can be addressed, using measurements at a small number of appropriately chosen baselines. The 
classification of Herbig Ae stars in two groups, based on the appearance of the spectral energy distribution (SED), 
has been attributed to a difference in disk geometry. Sources with a group I SED would have a flaring outer disk 
geometry, whereas the disk of group II sources is proposed to be flat (or "self-shadowed"). We show that this 
hypothesis can be tested using long-baseline interferometric measurements in the 10 fj,m atmospheric window. 
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1. Introduction 
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Herbig Ae/Be stars (HAEB Es, H erbig Il96(l for a more 
recent review see Natta et al. [loOOj are intermediate mass 
pre-main-sequence stars, surrounded by material which is 
left from the star formation process. A sub-group of mostly 
late B and A-F type HAEBE stars (hereafter HAEs) show 
little or no optical extinction, and usually have low mass 
ac cretio n rates as derived from radio analysis (Skinner et 
al. Il993j) and the lack of significant veiling in optical spec- 
tra. There is ample evidence that the circumstellar ma- 
terial responsible for the large infrared excesses of these 
stars is locate d in a circum stellar disk (e.g. Manni ngs fc 
Sargent Il997t Grad v et al. l200li Augereau et al. l200ll 
Eisner et al. l2003(l . Vink et al. l)2002l) show that the 
gaseous component has a disk-like geometry on scales of 
less than 0.1 AU. 

Whereas the presence of these circumstellar disks 
seems firmly established, the structure of the disks is a 
matter of debate. Kenyon & Hartmann l|l987f) developed 
"flaring" disk models for T-Tauri stars, in which H/R 
(the ratio of the disk surface height to the distance to the 
star) increases with increasing distance to the central star. 
The flaring disk model was refined by Chiang & Goldreich 
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l|l997l hereafter CG97) who introduced an optically thin 
surface layer responsible for the infrared emission fe atures 
generally seen in circ umste llar disks. Natta et al. (2001) 
and Dullemond et al. (2001) reconsider the CG97 model in 
the context of HAe stars, proposing that the innermost re- 
gion of the disk has an increased scale height: the "puffed- 
up inner rim". This configuration, which results from hy- 
drostatic equilibrium at the directly irradiated inner rim, 
naturally explains the near-infrared bump commonly ob- 
served in the Spectra l Ener gy Distribution (SED) of HAe 
systems (Natta et al. l200l|) . 

Meeus et al. l|200lh noted that based on the IR SED, 
HAEs can be divided into two main groups: "group I" 
sources that have a very strong, rising IR excess peaking 
around 60 /im, and "group II" sources displaying a more 
moderate IR excess, lacking the 60 /im bump. It was pro- 
posed that group I sources have a "flaring" geometry, al- 
lowing the disk to intercept and reprocess stellar radi- 
ation out to large stellocentric radii. In the outer disk 
of group II sources, on the other hand, H/R is approxi- 
mately constant, or decreasing with increasing distance to 
the star. The inner disk shields the outer disk from di- 
rect irradiation by the central star, hence the term "self- 
shadowed" disk. This substantially reduces the amount of 
radiation absorbed locally, leading to lower temperatures 
in the outer disk of a group II source. 

Recent 2 D mod eling by Dullemond (2002), Dullemond 
& Dominik ll2004 aj, henceforth DD04), has quantitatively 
confirmed that both flaring and self-shadowed disks are 
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natural solutions of the equation of vertical hydrostatic 
equilibrium in passive circumstellar disks (see also sec- 
tion 12. 4p . These models form the basis of the current 
study. There is ample evidence that group I and group II 
disks indeed have a flaring and s elf-sha dowed geome- 
try, resp ectively (e.g. Grady et al. 12004 Dullemond & 
Dominik r2004bl Leinert et al. l2004|) . However, as this is 
not an observational study, we will consistently refer to the 
models as flaring/self-shadowed, rather than group I /II 
(which is by definition an SED classification). 

With the advent of long (~ 10 2 m) baseline infrared in- 
terferometry using large apertertures, it has now become 
possible to observe HAe disks in the thermal infrared with 
a spatial resolution of order 10~ 2 arcsec. At the present, 
the number of baselines will be limited, and only interfer- 
ometric amplitudes (no phases) are available. True aper- 
ture synthetis imaging of disks is therefore not (yet) pos- 
sible. The interpretation of the measured visibility ampli- 
tudes, which contain information about the geometry of 
the disks, requires the use of disk models. 

In the near-infrared Herbig Ae/Be star disks have been 
observed wit h long-baseline interferometers since a num- 
ber of years jTuthill et all200ll iMillan-Gabet et "al~ll200lt 
lEisner et al1l2003ll2004 12005|) . Up to recently these mea- 
surements were evaluated using extremely simplified mod- 
els: Gaussian blobs, rings, ellipses etc. Such simple models 
made it possible to get a handle on the typical size and 
inclination of the emitting source, but did not go much fur- 
ther. For most sources the typical sizes were found to be 
in rough agreement with those predicted by the inner rim 
models, but at high accretion rates the ob servations de- 
viate from predictions. This is explained bv lAkeson et al.l 
( 200^) as due to the e mission from accretion inward of the 
inner dust rim, and bv lMonnier et all l)2005i) as due to the 
protection of dust by optically thick gas, allowing the dust 
to survive closer to the star. In these, and other, recent 
papers the modeling of the data already starts to go well 
beyond the simple ring/ellipse models, using actual multi- 
dimensional radiative transfer calculations in the case of 
Akeson et al., and detailed accr e tion d isk structure models 
in the case of ILachaume et all 1)200.^ . In particular with 
the new phase-closure capabilities in the near-infrared at 
the IOTA and VLTI/ AMBER interferometers such more 
advanced models are clearly of great use. 

With the mid-infrared interferometry capabilities of 
the MIDI instrument on the VLTI it is now possible to 
study the structure of the disk at slightly larger scales than 
the inner rim. This is the region in which the self-shadowed 
and flaring disks would most clearly be distinguishable. A 
first s et of measurements was published by iLeinert et alJ 
( 2004), and a first tentative correlation between the shape 
of the disk (flaring/self-shadowed) and the visibility was 
found. Mid-infrared interferometry also has the interesting 
capability of measureming mineralogical properties of the 
dust as a function of stellocentric radius. First measure- 
ments of this kind ( van Boekel et al1l2004a|) have revealed 
the strong radius-dependence of the crystallinity of dust, 
as predicted by theoretical models. In the present paper, 
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Fig. 1. A schematic representation of the geometries of 
the disk models studied in this work. The shadow cast by 
the inner rim is shaded light. 



however, we will be mostly concerned with the first aspect 
of mid-infrared interferometry: measuring the geometry of 
the disk. 

Bas ed on the 2D disk models of lDullemond fc Dominikl 
ll2004al) as well as the simpler models of 
lOhiang fc Goldreichl lll997ft we present calculations 
of the interferometric visibilities of HAe disks, to inves- 
tigate if it is possible to distinguish between the various 
disk geometries predicted by these theoretical models. 
Since the mid-infrared probes the structures at somewhat 
larger scale than the inner rim (from 1 AU out to about 
20 AU) , this wavelength regime is more suited to our aims 
than the near-infrared. The MIDI instrument is, so far, 
the only instrument capable of doing such measurements 
for Herbig Ae/Be stars, so in our analysis we focus on the 
typical baselines and properties of the VLTI. 



2. Modeling method 

The goal of this study is to predict and compare interfer- 
ometric visibilities of various disk models. Synthetic disk 
images are made using a ray tracing algorithm, where spe- 
cial care is taken to ensure all spatial scales in the disk are 
sufficiently resolved. Interferometric visibilities are calcu- 
lated by Fourier transforming the images. 
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2.1. The inner rim 

In the disk models considered here, the bright inner rim 
is treated in a highly simplified fashion: it is a sharp, 
"vertical wall". As a consequence, when such a model is 
viewed not pole-on, the flux from the near side of the in- 
ner rim is strongly suppressed since the hot, irradiated 
rim surface is occulted by the cooler parts immediately 
outside the rim surface. The far side of the bright inner 
rim is in full view and is responsible for essentially all of 
the near-infrared excess. In this configuration, the total 
amount of near-infrared emission depends strongly on the 
disk inclination, suggesting that the observed strength of 
the "3 micron bump" in the SED is a measure of the lat- 
ter. Observations of Herbig Ae stars, however, show that 
they all have rather similar near-infrare d exce sses, irre- 
spect ive of their inclination (Natta et al. 1200(1 Dominik 
et al. l2003|) . 

This indicates that the appearance of the inner rim is 
more smooth than the "vertical wall" used here. The pro- 
cesses that determine the shape of the inn er rim are cur- 
rently not yet understood. Isella & Natta 1)20051 ) recently 
showed that the dependence of the evaporation temper- 
ature on pressure naturally leads to a rounded-off inner 
rim. When such a disk is viewed at an inclination, both 
the near and the far side of the inner rim will be bright 
(although still the far side will be brighter). The bright 
inner rim will look like an inclined ring on the sky, rather 
than the "half ring" one obtains using the vertical wall 
model. 

Realistic radiative transfer modeling of a rounded-off 
inner rim introduces various numerical complications. To 
avoid these difficulties wc adopt the simplifie d vert ical rim 
structure used in Dullemond & Dominik (2004a). For a 
slightly off-polar inclination we artificially circularize the 
disk emission to circumvent the near-side/far-side asym- 
metry of the rim. In this way we mimic the rounded-off 
shape of the rim without having to confront the numer- 
ical complexities of radiative transfer in extreme optical 
depth rounded-off rims. While the spatial resolution of 
current 10 fim interferometers is just sufficient to measure 
the diameter of the inner rim, observations at higher res- 
olution are required to study details of the rim structure. 
We therefore believe that using this simplified approach is 
justified for our current purposes. To first approximation, 
inclination can be included by scaling the calculated spa- 
tial frequencies (or interferometric baselines) by a factor 
l/cos(#) along the minor axis of an inclined disk, where 
6 is the inclination of the disk. At high inclinations, this 
approximation brakes down. 

2.2. Interferometry 

Rather than images, an interferometer produces an inter- 
ference signal called the interferometric "visibility" (V), 
which is the spatial coherence function of the intensity 
distribution of the source. The visibility is related to the 
intensity distribution through the van-Cittert-Zernike the- 
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Fig. 2. The spectral energy distributions of the DD04 flar- 
ing (full curve), DD04 self-shadowed (dashed curve) and 
CG97 (dotted curve) models. The full grey curve repre- 
sents the DD04 flaring model, where the silicate reso- 
nances between 8 and 25 /im have artificially been removed 
from the opacity table (see also section 13 . 3 . 21 and Fig.EJ. 

orem, which states that the visibility is the Fourier trans- 
form of the intensity distribution of the source. For an 
introduc tion t o long-baseline interferometry, we refer to 
Lawson l|2000h . 

2.3. Spectrally resolved visibilities 

Traditionally, visibility curves are represented as a func- 
tion of the projected baseline B 1 , at a specific wavelength 
(a V(B) curve). Such a curve represents a number of visi- 
bility measurements at different projected baselines, which 
usually requires the use of multiple telescope pairs and/or 
moving telescopes. 

A new possibility in the 10 fj,m region is the use of spec- 
tral dispersion with large wavelength coverage. When us- 
ing an instrument that has this capability, one can obtain a 
whole "visibility curve" in one single measurement. Unlike 
the common V(B) curve, the V(A) curve obtained this 
way holds many visibility values at only one baseline. The 
spatial resolution of the observation (ss B/X) changes by 
almost a factor of 2 between 7.5 and 14 /im. Detailed mod- 
eling is required to interpret V(A) curves. Most HAe stars 
show a prominent emission band between 8 and 12 pm, 
due to silicate dust. The shape of this emission band varies 
strongly, depending on chemical composition, particle size 
and lattice structure of the silicate grains. When simulat- 
ing interferometric visibilities using disk models, one finds 
that the detailed shape of the visibility curve depends on 
the opacities used, i.e. on the dust properties. These are 

1 The spatial frequency of the observation is: k = B/X, where 
A is the observing wavelength. The units of k are cycles/radian 
if A and B have the same units 
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Fig. 3. Radial intensity profiles of the DD04 flaring, DD04 
self-shadowed and CG97 models (upper panel). The lower 
panel shows the normalized cumulative flux distribution 
of the models. The three regions of the disk that we dis- 
tinguish ( "bright inner rim" , "intermediate shadowed re- 
gion" , and "outer disk" ) are indicated for the DD04 flaring 
model (in the DD04 self-shadowed model, the intermedi- 
ate shadowed region extends somewhat further outward). 



different from star to star, and vary within a disk as a 
function of distance to the star l(van Boekel et all ferKMa - ! . 
The visibilities measured in the silicate emission feature 
are a mixture of disk structure and mineralogy. Therefore, 
in order to deduce information on the disk structure, it is 
preferable to use measurements at wavelengths outside the 
silicate feature, which in practice means between about 12 
and 13.5 /im. At 8 /im, it is also possible to sample the con- 
tinuum emission, this is however more difficult since here 
the atmospheric transmission is rather poor. 



we use the following stellar parameters: M* = 2.5 M©, 
i?„ = 2R© and T» = 10 000 K, which amounts to a stellar 
luminosity of L* = 36 Lq. All disk models in this work 
have a disk mass of 0.1 M©, a gas-to-dust ratio of 100, a 
surface density distribution of E oc i? -15 , and an outer 
radius of 200 AU. The DD04 models have an inner disk 
radius which is calculated self-consistently assuming an 
optically thick inner rim. The location of the inner rim is 
set by the dust evaporation temperature, which is about 
1500 K for silicate dust. The CG97 model has an inner 
radius of 0.21 AU, which corresponds to the radius where 
the black-body temperature is 1500 K. For the dust opaci- 
ties we use a simple m odel co nsisting only of small silicate 
grains (Laor & Draine ll998h . 

If the disk is optically thick en ough, the disk has a 
flaring shape (DD04 and Dullemond l2002l) . When the op- 
tical depth is decreased, a flaring disk can turn into a 
self-shadowed disk and the SED changes from a group I 
to a group II shape. The flaring and self shadowed models 
shown here are the BL1 and the BL4 model from DD04. 
Both disks have a mass of 0.1 M©. In the BL1 model all 
the dust mass is in 0.1 /im silicate grains. In the BL4 model 
99.9% of the mass has been converted into 2 mm size grains 
located in the midplane, while only 0.1% remains in small 
0.1 /im grains, thus strongly lowering the opacity of the 
disk. We stress however, that the flaring vs. self-shadowed 
behaviour of the disk depends on high vs. lower optical 
depth, and that dust coagulation is a possible mechanism 
to achieve lower optical depths. 

In the proposed scheme, the outer disk of a self- 
shadowed source is shielded from direct stellar irradiation 
by its own inner disk. However, the outer disk receives 
near-IR radiation emitted by the hot innermost disk re- 
gions, and optical/UV radiation which is scattered by the 
diffuse inner disk atmosphere. Therefore, the temperature 
and scale height in a self-shadowed disk are still signifi- 
cantly larger than zero. Note that in a flaring disk, there is 
also a region just outward of the inner rim that is shielded 
from direct stellar radiation. Contrary to a self-shadowed 
disk, a flaring disk emerges from the shadow cast by the 
puffed up inner rim, at distances of a few AU from the 
star. 



2.4. The DD04 disk model 

The disk models used in this work are described in DD04. 
These are 2D axisymmetric models in which the gas and 
dust density and temperature are given as a function of 
radius R and polar coordinate 0. The disk is assumed to 
be heated only by irradiation by the central star. A 2D 
continuum radiative transfer code is used to compute the 
entire temperature structure of the disk. The vertical den- 
sity structure, for a given radial surface density distribu- 
tion E(i2), is computed by demanding vertical hydrostatic 
equilibrium. In this way the disk has a self-consistent tem- 
perature and density structure, from which images and 
SEDs can be computed using a ray-tracer. For this work 



3. Results 

The geometry of the disk models investigated in this work 
is schematically represented in Fig. ^ The outer disk ge- 
ometries in the CG97 and DD04 flaring models are very 
similar, however the inner disk region is different. The 
DD04 flaring model has a puffed up inner rim, and an 
intermediate shadowed region (shaded light). The DD04 
self-shadowed model has an inner disk structure that is 
very similar to the DD04 flaring model. The outer disk 
in the self-shadowed model never rises above the shadow 
cast by the puffed up inner rim, but is irradiated by the 
hot inner disk regions. 
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3.1. Spectral energy distributions 

The emerging spectral energy distributions of the disk 
models are shown in Fig. [2] The infrared excess emis- 
sion at far infrared (FIR) wavelengths (^60-100 /im) is 
clearly stronger in the flaring disk models than in the self- 
shadowed model. The FIR excess is significantly stronger 
in the DD04 flaring model than in the CG97 (also flaring) 
model. This can be traced to the simplifications made in 
the CG97/DDN01 models which do not take properly into 
account various 3D radiative transfer effects. In particular, 
these models do not account for the moderate 'boosting' 
of radiation toward the polar axis to compensate for the 
occultation in equatorial directions by the disk's own flar- 
ing outer regions. In the DD04 models these effects are 
consistently taken into account by virtue of the full multi- 
dimensional radiative transfer treatment used in those 
models. The near-infrared excess around 2-3 /im which is 
very prominent in the DD04 model SEDs, is much less 
pronounced in the CG97 model. In the 10 /im region, the 
CG97 and DD04 flaring model SEDs are virtually iden- 
tical in spectral shape, though the DD04 flaring model 
has a somewhat higher absolute flux level. The DD04 self- 
shadowed disk is fainter than the flaring models, and has 
a bluer continuum slope in the 10 /im region. 

3.2. Radial intensity profiles 

Fig. El shows the radial intensity profiles at 12.6 /zm of a 
DD04 flaring and self-shadowed disk model, and the CG97 
model. Both the DD04 flaring and self-shadowed model 
essentially exhibit three regimes: 

1) the "bright puffed up inner rim", that causes a ringlike 
emission, contributing mainly between 0.5 and 0.8 AU 
from the star. 

2) a region just behind the inner rim (as seen from the 
star), where the dust temperatures are much lower 
than in the inner rim. From this "intermediate shad- 
owed region", relatively little radiation emerges (as 
can be seen in the cumulative flux distributions of the 
DD04 models in the lower panel of Fig. [31 which are 
nearly constant in this region). 

3) the "outer disk region" , whose main flux contribution 
arises between 3 and 20 AU from the central star. In 
a flaring disk model, this outer disk region is directly 
irradiated by the central star. 

The "bright puffed up inner rim" , "intermediate shad- 
owed region" and "outer disk region" are of course just 
different parts of the same physical structure, and the 
distinction made here serves merely to help the reader 
develop a qualitative understanding of how such geome- 
tries translate into interferometric visibilities. In Fig.|31we 
have indicated the three regions discussed above. For the 
DD04 flaring model, the outer radius of the intermediate 
shadowed region can be well defined to be between 2 and 
3 AU, where the slope of the cumulative flux distribution 
clearly increases. In the DD04 self-shadowed model, this 



radius is less clearly defined but evidently somewhat larger 
than in the DD04 flaring model. The contribution of the 
bright inner rim emission to the total system flux depends 
strongly on wavelength (for the self-shadowed model, the 
inner rim contributes more than 90% to the total flux at 
6 /im, about 60% at 8 /im, about 35% at 13 /im and less 
than 5% at 30 /im) , and is always higher in a self-shadowed 
model than in a flaring model harboring the same central 
star. 

The CG97 model has, per definition, a flaring disk 
structure. Contrary to the DD04 flaring model however, 
it does not have a bright puffed up inner rim, and conse- 
quently it lacks an intermediate shadowed region. 

3.3. Visibility curves 

3.3.1. Visibility as a function of baseline 

Fig. 01 shows the predicted "classical" visibility curves 
V(B) of the considered disk models, at several wave- 
lengths. From top to bottom we show the predictions for 
the CG97 model, the DD04 flaring model and the DD04 
self-shadowed model. To develop an intuitive understand- 
ing of how the characteristics of the emerging intensity 
distributions of the various models are reflected in their 
visibility curves, we will discuss the curves in Fig. 01 In this 
example, the star is put at a distance of 150 pc, typical for 
nearby Herbig stars. 

At a baseline of m all sources are of course unresolved 
and have a visibility of 1. The CG97 model shows a steady 
drop in visibility as the baseline is increased. The slope 
of the visibility curve changes gradually, reflecting that 
the radial intensity profile shown in Fig.El has no strong 
substructure. In this sense, the CG97 model is "scaleless" . 

For the DD04 models, this is different. These models 
have essentially two scales: the bright inner rim, which 
emits between 0.5 and 0.8 AU, and the outer disk, which 
emits most of its flux between about 3 and 20 AU. In be- 
tween lies the intermediate shadowed region, whence little 
flux emerges. This general picture is reflected in the visi- 
bility curves. Starting at m, and increasing the baseline, 
we observe a steady drop in visibility as the outer disk 
gets more and more resolved. Note that at a baseline of 
10 m the visibility is already significantly lower than 1, 
predicting that the largest modern day telescopes might 
marginally resolve the outer disks in such objects at 10 /zm. 
For the H Ae star HD 100546 this ha s indee d been observed 
(Liu et al. I2003L van Boekel et al.|2Q04b). At a baseline 
of about 30 m, the outer disk is mostly resolved while the 
inner rim is still essentially unresolved. Therefore, the vis- 
ibility curves flatten at this point. The visibility level at 
this baseline (about 15% for the flaring model and 40% 
for the self-shadowed disk, at 9.8 /im) indicates the frac- 
tion of the total system flux that is emitted by the bright 
inner rim. At longer baselines, the bright inner rim itself 
becomes resolved by the interferometer, and the visibility 
gradually goes to its first null. Since the spatial resolu- 
tion of the interferometer scales inversely with wavelength, 
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Fig. 4. Simulated visibility curves V(B) of a CG97, a 
DD04 flaring, and a DD04 self-shadowed disk model (from 
top to bottom, for V(A) curves of these models see Fig.EJ). 
We show model visibilities at three different wavelengths. 
x\ and #2 denote optimum baselines to distinguish be- 
tween the various models (see section . 



Fig. 5. Visibility curves V(A) of a CG97, a DD04 flaring, 
and a DD04 self-shadowed disk model (for V(B) curves 
of these models see Fig.0J. We show predicted behaviour 
of the visibility as a function of wavelength V(X), for sev- 
eral different baselines. The wavelength regions that are 
inaccessible from the ground are shaded grey. 



whereas the apparent diameter of the inner rim hardly de- 
pends on the wavelength, zero visibility is reached first at 
the shortest wavelengths, and at longer baselines for the 
longer wavelengths. 



3.3.2. Visibility as a function of wavelength 

Fig.EOshows V(X) curves for the CG97, DD04 flaring and 
DD04 self-shadowed model, at a number of different base- 
lines. Each of these curves represents a single, dispersed, 
visibility measurement (the spectral region inaccessible 
from the ground is shaded grey). Each curve in Fig.[3can 
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be regarded as a cut through Fig.0|at a specific baseline, 
with a much denser wavelength sampling. 

The overall trend for all curves is to show the highest 
visibilities at 8 /im, and lower visibilities at 13 /im. This 
is because the apparent size of the disks increases with 
wavelength more rapidly than the intcrferometric resolu- 
tion decreases. There is generally a sharp decrease in vis- 
ibility between 8 and 10 ^m. There are two reasons for 
this. For the models with an inner rim one reason is that 
the emission from this rim dominates the spectrum below 
about 8 ^m. The emission at 10 /jm originates from more 
extended regions of the disk, resulting naturally in a lower 
visibility than the 8 /im emission. A second reason for the 
decline of the visibility between 8 and 10 /im - and for the 
gentle rise in visibility toward 13 /jm in some models - is 
that the flux in the 10 /im silicate feature originates pre- 
dominantly from the warm surface layers of the disk, while 
the flux outside the feature comes from the cooler regions 
below. The warm dust in the surface layer can radiate in 
the mid-infrared out to larger radii than the cooler dust in 
the disk interior. In other words: in the warm surface lay- 
ers the Wien exponential cut-off in the mid-infrared takes 
effect at larger radii than in the disk interior. This expla- 
nation also holds for the CG97 models, which do not have 
an inner rim. 

In FigEJwe demonstrate the importance of the silicate 
resonance for the simulated visibility curves. We show the 
visibilities of the DD04 flaring model (full curves, see also 
the middle panel of Fig.[5j). To calculate visibility curves 
of a model without silicate resonances, we removed the 10 
and 20 micron silicate features from the opacity table prior 
to the ray tracing (dotted curves). The large influence of 
the opacity of the material on the resulting visibilities is 
evident. The SED of the model without silicate resonances 
is indicated by the full grey curve in Fig.|3 

The interpretation of the curves in Fig.[5] in terms 
of disk geometry is not straightforward, for several rea- 
sons. First, the emerging intensity distribution of the disk 
changes with wavelength. Second, the spatial resolution of 
the interferometer decreases by almost a factor of two be- 
tween the short and long wavelength edges of the 10 fim at- 
mospheric window (N-band). Third, how the disk intensity 
distribution (and thus the visibilities) changes with wave- 
length depends on the opacity of the dust, i.e. on mineral- 
ogy (Fig. El). When we measure a single V(X) curve, what 
we see therefore is a mixture of disk geometry changing 
with wavelength, instrumental resolution changing with 
wavelength, and the mineralogy of the source. 

As the mineralogy changes from star to star, it is diffi- 
cult to obtain general diagnostics for disk geometry from 
a V(X) curve. The detailed interpretation of such mea- 
surements requires a model of each individual star, where 
both the spectrum (mineralogy) and the disk structure 
are fitted simultaneously. As a general diagnostic for disk 
structure it is thus preferable to measure visibilities in the 
continuum, where the visibilities do not depend strongly 
on mineralogy. We will come back to this in section l3~4l 




10 12 
A [/*m] 

Fig. 6. The effect of the increased dust opacity in the 
silicate feature on the visibilities. The full lines show the 
visibility curves of the DD04 flaring model, the dashed 
lines show the visibility curves of the same model, where 
the silicate resonances have been removed in the opacities 
prior to the ray tracing. 



However, the spectral capabilities of the interferometer 
develop their full strength, when one uses observations at 
several baselines in order to reconstruct an "image" of the 
disk 2 . With spectrally dispersed visibilities, the spectrum 
of the disk is then known immediately at all positions in 
the disk. It is then possible to study the mineralogy, size 
distribution and chemical composition of dust grains in 
the disk surface layer as a function of distance to the cen- 
tral star, providing crucial information about processing 
and radial mixing in disks. 

If the intensity distribution in the disk is strongly cen- 
trally peaked like the models discussed in the present pa- 
per, the correlated flux obtained at a single, long base- 
line 3 can be directly interpreted as the spectrum of the 
innermost regions of the disk (in the correlated spectrum 
obtained with only one measurement, there is still an un- 
known spatial term mixed in, that typically introduces a 
slope in the spectrum. This however has little influence 
on the derived mineralogy). The outer disk spectrum can 
then be obtained as a difference between the integrated 
disk spectrum and the inner disk spectrum. Applying this 
method to the first spectrally resolved full N-band visibil- 



2 In practice this will be easiest for disks that are seen not too 
far from pole on, so their image has a high degree of azimuthal 
symmetry. This limits the number of baselines needed, and 
requires the measurement of visibility amplitude only, since all 
phases will be approximately 

3 The visibility is the ratio of the correlated flux and the total 
flux (F cor /Ftot). The correlated flux, or correlated spectrum in 
the case of a spectrally dispersed measurement, is the quantity 
an interferometer measures. 
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Fig. 7. A diagnostic diagram used to distinguish between 
models with (DD04) and without (CG97) a bright inner 
rim, and between models with a flaring (DD04 flaring) and 
self-shadowed (DD04 self-sh) outer disk geometry. On the 
horizontal axis we plot the slope of the visibility curves 
between two appropriately chosen "normalized" baseline 
lengths X\ and X2 (x* = B i /d 1 where Bi is the baseline in 
m and d is the distance to the star in pc, see section I3~4l 
for how X\ and x-i are best chosen). On the vertical axis 
we plot the predicted visibility at the longest baseline. 
The CG97 models have a much steeper slope than the 
DD04 models. The DD04 self-shadowed model has a much 
higher visibility than the DD04 flaring model. In the upper 
left corner we have indicated the uncertainties due to the 
limited precision of the visibility measurements, where we 
have assumed a 1% accuracy in visibility. In grey symbols 
we have plotted where in the diagram the various models 
end up if we artificially remove the silicate feature (as we 
did in Fig. [fjj). This can be regarded as an upper limit for 
the uncertainty due to mineralogy. 

ity measurements of HA e stars, it was demonstrated by 
Ivan Boekel et alJ l|2004a[) that the mineralogy in the disk 
can vary strongly with distance to the star. 

3.4. Distinguishing between the various models 

The goal of this study is to show how interferometric mea- 
surements can be used to distinguish between the various 
disk models. Clearly, the curves in Fig.QJare different for 
the different models. However, one will typically not have 
continuous measured visibility curves at hand, but rather 
have samples at a few different baselines. Here we show 
that it is, at least in principle, possible to distinguish both 
between the CG97 and DD04 models on one hand, and be- 
tween the DD04 flaring and DD04 self-shadowed on the 
other, using measurements at only two appropriately cho- 
sen baselines. 



The distinction between the CG97 and DD04 mod- 
els is based on the absence of an intermediate shadowed 
region in the former. At spatial scales corresponding to 
the intermediate shadowed region in the DD04 models, 
little flux emerges. Therefore, the visibility curves are rel- 
atively flat at the baselines corresponding to these spa- 
tial scales, they show a "plateau" (very prominent in the 
DD04 flaring model visibility at 12.6 in Fig.0J. The 
CG97 model does not have such specific spatial scales 
with much reduced emergent intensity, and therefore lacks 
the plateau in the visibility curve. The visibility curves of 
CG97 and DD04 models thus have a different slope at 
baselines corresponding to the scale of the intermediate 
shadowed region. Once this slope difference has been de- 
tected, DD04 flaring and DD04 self-shadowded disks can 
be distinguished by the relative contribution of the bright 
inner rim to the total system flux, which is much higher 
for a self-shadowed model. Note that the bright inner rim 
itself is virtually identical in both models, but the outer 
disk is much brighter in the flaring disk than in the self- 
shadowed case. Therefore, the predicted visibilities at our 
selected baselines are much lower for the flaring model. We 
recall that deducing properties about the disk structure is 
best done outside the silicate emission feature, which in 
practice favours the region between 12 and 13 /im. 

A measurement at a specific baseline samples the cor- 
responding angular scale, and the physical scale (in AU) 
associated with this baseline therefore depends linearly on 
the distance to the star. It is therefore convenient to in- 
troduce the "normalized baseline" 



where B is the baseline in meter and D is the distance to 
the star in parsec. Consider the emission of a small part of 
the disk, arising between R and R + dR from the central 
star, i.e. a ring with angular diameter 9 — R/d (where d 
is the distance to the star) . The visibility curve of such an 
annulus of emission is the zeroth order Bessel function: 

R 

V an (B) = J ( — ) (2) 

where B is the interferometric baseline, and A is the wave- 
length of observation. The visibility reaches the first null 
at a baseline of about 

B « 158-^ [m] (3) 

"mas 

where for convenience the wavelength and annulus angular 
diameter have been expressed in fim and milli-arcseconds, 
respectively. Our goal is to detect the effect of the inter- 
mediate shadowed region, i.e. the very low flux contribu- 
tion from annuli between about 0.8 and 3 AU. To estimate 
which baseline is most sensitive to emission from an an- 
nulus with diameter 0, let us take the baseline where the 
visibility has half its maximum value, U an =0.5: 

#0.5 « 100^ [m] (4) 
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For an annulus of 2 AU radius, this corresponds to a nor- 
malized baseline of x = A Mm /40. Let us take the region 
between 1 and 2 AU from the star as characteristic for the 
intermediate shadowed region. 

We find that in order to most clearly separate the spe- 
cific DD04 models (with an intermediate shadowed region) 
used in this work from the CG97 models (without an in- 
termediate shadowed region), the best choice for x is: 



x l = 0.038 

X2 = 2xi 



(lO Mm) (lq) 



(5) 



To allow the above formula to be applied to stars of 
different luminosity we have applied a simple scaling of 
x with the square root of the luminosity. The visibility 
is now sampled at the baselines corresponding to x\ and 
X2 (for our 36 star at 150 pc and at a wavelength of 
12.6 /im these are 43 and 86 m, respectively, see also Fig0J), 
yielding visibilities V\ and V2 ■ In Fig.0 we plot the on 
the vertical axis the visibility measured on the longest 
of the two baselines (V2). On the horizontal axis we plot 
the deduced slope of the visibility curve between X\ and 
X%. In addition to the CG97 model with an inner radius 
of 0.21 AU we plot CG97 models with inner radii of 0.02 
and 0.5 AU, to illustrate the behaviour with varying inner 
radius. From this figure we see that: 

1) The CG97 models have a much steeper slope than the 
DD04 models in this baseline regime. 

2) The DD04 flaring model has a much lower visibility at 
the long baseline than the DD04 self-shadowed model 
(and on the short baseline as well). 

We can conclude therefore, that it is possible to discrimi- 
nate between the various models, with a very limited num- 
ber of measurements. In the upper left corner of Fig.0 
we indicate the uncertainty in the measured location of 
a source in this diagram due to limited precision of the 
visibility measurments. We have here assumed a visibility 
accuracy of 1 %, which is what MIDI aims to achieve. The 
difference in the position of the models with and without 
a silicate feature (plotted with black and grey symbols, 
respectively) can be regarded as the extreme case for the 
uncertainty arising from mineralogy, since for each indi- 
vidual source, we know what the silicate feature looks like. 

This analysis has been done in an idealized world 
where we have both assumed that the disks are pole-on, 
and that our models are a good representation of the true 
disk geometry. The optimum choice for X\ and xi depends 
on the geometry of the disk and may therefore in reality 
be somewhat different than the values given in our recipe 
(equation El . If the outer disk is smaller than we predict, 
the value of x\ should be increased. If the bright inner 
rim is located at radii somewhat larger than predicted, the 
value of xi should be decreased. There is however evidence 
that the bright inner rim is located at radii somewhat 
smaller than predicted in our models (Eisner et al. 2003), 
and therefore this is not likely a reason for concern. In 



practice, measurements at more than two baselines (~5) 
are probably needed to unambiguously establish the na- 
ture of the sources. 

4. Conclusions 

We have presented model calculations of interferometric 
visibilities of circumstellar disks around Herbig Ae stars. 
We compare predictions for disks with (DD04) and with- 
out a bright inner rim (CG97). We show that it is possi- 
ble to distinguish between both model possibilities using 
a small number of interferometric measurements in the 
10 /im atmospheric window. Such measurements also allow 
to distinguish between flaring and self-shadowed disk mod- 
els. This allows testing of the hypothesis that group I and 
group II sources correspond to flaring and self-shadowed 
disks, as is suggested by their spectral energy distribu- 
tions. 

Acknowledgements. R. Lachaume and L. B. F. M. Waters are 
gratefully acknowleged for constructive comments on the 
manuscript. 

References 

Akeson, R. L., Walker, C. EL, Wood, K., et al. 2005, ApJ, 
622, 440 

Augereau, J. C., Lagrange, A. M., Mouillet, D., & Menard, 

F. 2001, A&A, 365, 78 
Chiang, E. I. & Goldreich, P. 1997, ApJ, 490, 368 
Dominik, C, Dullemond, C. P., Waters, L. B. F. M., & 

Walch, S. 2003, A&A, 398, 607 
Dullemond, C. P. 2002, A&A, 395, 853 
Dullemond, C. P. & Dominik, C. 2004a, A&A, 417, 159 
— . 2004b, A&A, 421, 1075 

Dullemond, C. P., Dominik, C, & Natta, A. 2001, ApJ, 
560, 957 

Eisner, J. A., Hillenbrand, L. A., White, R. J., Akeson, 

R. L., & Sargent, A. I. 2005, ApJ, 623, 952 
Eisner, J. A., Lane, B. F., Akeson, R. L., Hillenbrand, 

L. A., & Sargent, A. I. 2003, ApJ, 588, 360 
Eisner, J. A., Lane, B. F., Hillenbrand, L. A., Akeson, 

R. L., & Sargent, A. I. 2004, ApJ, 613, 1049 
Grady, C. A., Polomski, E. F., Henning, T., et al. 2001, 

AJ, 122, 3396 

Grady, C. A., Woodgate, B., Torres, C. A. O., et al. 2004, 

ApJ, 608, 809 
Herbig, G. H. 1960, ApJS, 4, 337 
Isella, A. & Natta, A. 2005, in press 
Kenyon, S. J. & Hartmann, L. 1987, ApJ, 323, 714 
Lachaume, R., Malbet, F., & Monin, J.-L. 2003, A&A, 

400, 185 

Laor, A. & Draine, B. T. 1993, ApJ, 402, 441 

Lawson, P. R., ed. 2000, Principles of Long Baseline Stellar 

Interferometry 
Leinert, C, van Boekel, R., Waters, L. B. F. M., et al. 

2004, A&A, 423, 537 
Liu, W. M., Hinz, P. M., Meyer, M. R., et al. 2003, ApJ, 

598, Llll 



10 



van Boekel et al.: Herbig Ae star disks with a 10 /jm interferometer 



Mannings, V. & Sargent, A. I. 1997, ApJ, 490, 792 
Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, 

A&A, 365, 476 
Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, 

ApJ, 546, 358 

Monnicr, J. D., Millan-Gabet, R., Billmcicr, R., et al. 

2005, ApJ, 624, 832 
Natta, A., Grinin, V., & Mannings, V. 2000, Protostars 

and Planets IV, 559 
Natta, A., Prusti, T., Neri, R., et al. 2001, A&A, 371, 186 
Skinner, S. L., Brown, A., & Stewart, R. T. 1993, ApJS, 

87, 217 

Tuthill, P. G., Monnicr, J. D., & Danchi, W. C. 2001, 

Nature, 409, 1012 
van Boekel, R., Min, M., Leinert, C, et al. 2004a, Nature, 

432, 479 

van Boekel, R., Waters, L. B. F. M., Dominik, C, et al. 

2004b, A&A, 418, 177 
Vink, J. S., Drew, J. E., Harries, T. J., & Oudmaijcr, R. D. 

2002, MNRAS, 337, 356 



